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ABSTRACT 

We present the first high resolution X-ray image of the dwarf elliptical galaxy NGC 3226. The data 
were obtained during an observation of the nearby Seyfert Galaxy NGC 3227 using the Chandra X-ray 
Observatory. We detect a point X-ray source spatially consistent with the optical nucleus of NGC 3226 
and a recently-detected, compact, fiat-spectrum, radio source. The X-ray spectrum can be measured 
up to ~10 keV and is consistent with a power law with a photon index 1.7 < T < 2.2, or thermal 
bremmstrahlung emission with 4 < kT < 10 keV. In both cases the luminosity in the 2-10 keV band 
~ 10 40 /i 

75 er g s 1 - We find marginal evidence that the nucleus varies within the observation. These 
characteristics support evidence from other wavebands that NGC 3226 harbors a low-luminosity, active 
nucleus. We also comment on two previously-unknown, fainter X-ray sources < 15 arcsec from the 
nucleus of NGC 3226. Their proximity to the nucleus (with projected distances < I.3/I75 kpc) suggests 
both are within NGC 3226, and thus have luminosities (~fewx 10 38 -fewx 10 39 erg s _1 ) consistent with 
black- hole binary systems. 

Subject headings: galaxies: dwarf - galaxies: active - galaxies: individual (NGC 3226) - galaxies: 
nuclei - X-rays: galaxies - X-rays: binaries 



1. INTRODUCTION 

NGC 3226 is a dwarf elliptical galaxy (z = 0.00441 ± 
0.00006; e.g. de Vaucouleurs et al. 1991) interacting with 
the nearby SAB(s)a galaxy NGC 3227 9 . The center of 
NGC 3226 hosts a low-ionization nuclear emission-line re- 
gion (LINER) of optical spectroscopic type 1.9 (Ho, Fil- 
lipcnko & Sargent, 1997a), along with large quantities of 
dust in filamentary structures (e.g. Rest et al. 2001). 
A compact, flat-spectrum radio source has recently been 
detected (Nagar et al. 2000; Falcke et al. 2000) at a loca- 
tion consistent with the optical nucleus. This, along with 
the claim of a broad component to the Ha emission line 
(FWHM ~ 2 x 10 3 km s" 1 ; Ho et al. 1997b), suggests the 
galaxy hosts an active galactic nucleus (AGN) rather than 
a region of aging starburst activity. The luminosity of the 
broad Ha line in NGC 3226 (L(Ha) ~ 10 39 erg s" 1 ) is 
an order of magnitude below the dividing line between so- 
called low-luminosity AGN (LLAGN, or "dwarf" Seyfert 
galaxies) and the more powerful ( "normal" ) Seyfert 1 and 
giant elliptical galaxies. 

X-ray emission from the nuclear regions of NGC 3226 
has been detected using both the .RO&4TPSPC (Komossa 
& Fink 1997; Sansom, Hibbard & Schweizer 2000) and 



HRI (Roberts & Warwick 2000) with an implied luminos- 
ity L(0.1-2 keV) ~ 10 40 erg s _1 . However these data sets 
had relatively poor spatial resolution, and were unable to 
distinguish either spatially or spectroscopically between 
an (AGN-like) power law X-ray spectrum and a collision- 
ally ionized plasma due to hot gas (e.g. from a starburst 
region) . 

In this Letter we present results for NGC 3226 obtained 
using the Chandra X-ray Observatory (CXO: e.g. Weis- 
skopf, O'Dell, van Speybroeck 1996). Following a brief 
description of the observation (§2), in §3 we report the 
detection of three X-ray sources within the inner ~ 2 kpc 
of the nucleus, with brightest source identified as the nu- 
cleus of NGC 3226. In §4 we report the temporal and 
spectral characteristics of the nucleus, and discuss the ev- 
idence that NGC 3226 harbors a LLAGN. The limited 
constraints that can be obtained from the two off-nuclear 
sources are discussed in §5. We present our conclusions in 
§6- 

2. THE OBSERVATION 



The data from NGC 3226 reported here were obtained dur- 
ing an observation using the CXO in 1999 Dec. The High 
Energy Transmission Grating Spectrometer (HETGS: e.g. 
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Fig. 1. — The O th -order HETGS image of the nuclear region of NGC 3226 in the 0.8-6 keV band. For clarity the image has been smoothed 
by a Gaussian with a = 0.492 arcsec. The DSS image is overlaid as the blue, linearly-spaced contours. The separation between the grid lines 
(10 — 2 degrees) ~ kpc. The offsets between the position of the brightest X-ray source and the optical nucleus of NGC 3226 as tabulated 

by Cotton, Condon & Arbizzani (1999) are (ARA, Adec) = (—1.2 arcsec, +0.2 arcsec). Currently there are no optical identifications of the 
other two X-ray sources (marked X-2 & X-3). 



Markert et al. 1994) was employed with the Advanced 
CCD Imaging Spectrometer (ACIS: e.g. Nousek et al 
1998) in the focal plane (with a temperature of —110° C). 
The HETGS consists of two sets of gratings (hereafter 
the medium- and high-energy grating {MEG & HEG} 
"arms"), each with a different period and intercepting X- 
rays from (two of the four) different shells of the X-ray 
telescope. The rulings of the MEG and HEG differ by 
~10 degrees preventing any spatial overlap in the dispersed 
spectra of an on-axis, point source. The ACIS-S array was 
placed at the optical axis, consisting of 6 CCD chips 10 
orientated to record the undispersed (0 th -order) photons 
from the target and nearby sources (on the ACIS-S3 CCD), 
along with both the positive and negative orders of the 
spectrum dispersed by both grating arms of the HETGS. 

The observation was continuous, with a total duration 
of 49.9 ks. Following screening and removing the "streaks" 
from the s4 chip using destreak (v 1.3; Houck 2000), the 
resultant exposure time was 46.3 ks. All the data analysis 
presented here was performed using the Chandra CALDB 
(v2.0), and the CIA0 (v2.0.2) and HEAsoft (v5.0.2) soft- 
ware packages. Further details on the instrumental param- 
eters and the data reduction can be found in George et al. 
(2001), along with the scientific results for NGC 3227. 



3. SPATIAL ANALYSIS 



In Fig. 1 we show an image of NGC 3226 in the 0.6- 
8 keV band. For simplicity, here we ignore any pecu- 
liar motion associated with the NGC 3226/3227 system 
or other members of its local group, and simply use the 
heliocentric velocity of NGC 3226 to derive a distance 
to 17.6/ifg 1 Mpc (where h 75 = H /75 km s" 1 Mpc -1 ). 
Thus 1 arcsec = 85hj^ pc. Three X-ray sources are evi- 
dent. We find the centroid of the brightest source to be 
at RA=10h 23m 27.00s, dec=+19d 53m 54.7s (J2000). 
This is consistent the position of the optical nucleus of 
NGC 3226 (RA=10h 23m 27.08s, dec=+19d 53m 54.5s; 
Cotton, Condon & Arbizzani 1999) within the cur- 
rent uncertainties in the aspect reconstruction of CXO 
data (<1 arcsec; Aldcroft, p.comm), and in excellent 
agreement with the compact, flat-spectrum radio source 
(RA=10h 23m 27.01s, dec=+19d 53m 54.5s; Falcke et al. 
2000). Thus we identify this X-ray emission as the HETGS 
th -order image of the inner - 150/ifg 1 pc of NGC 3226. 
From an analysis using a Marr wavelet function provided 
by the wavdetect task in CIA0 we estimate a net source 
count rate of 491 ± 22 count s" 1 . 

To the best of our knowledge, nothing has been re- 
ported in the literature concerning the other two sources, 
CXO J102334. 1+195347 and CXO J102326. 7+195407 (la- 
belled X-2 & X-3 respectively in Fig. 1). Both are much 
weaker than the nucleus (with estimated source count rates 
of 20 ± 5 count s _1 and 18 ± 4 count s _1 respectively). 



The ACIS-S0 chip was not in operation during the observation due to instrumental concerns. 
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Fig. 2. — a) A light curve of the nucleus of NGC 3226 in the 0.8-6 keV band using 2 X 10 3 s bins. Note that the y-axis is logarithmic, with a 
dynamic range 0.33-3.0 times the mean count rate (0.0099 ct s — 1 ; dashed line). The predicted background rate is ~ 3 X 10 -5 count s -1 and 
constant, b) Observed distribution (red) of number of bins (each of 10 s) containing n counts, and a Poisson distribution (blue) derived from 
the mean count rate. The dashed lines show the 95% confidence limits for the Poisson distribution. The bins marked with an "X" highlight 
more bins are observed outside these limits than predicted. 



4. THE NUCLEUS OF NGC 3226 

For the (temporal and spectral analysis) presented below 
we employed a circular extraction cell of radius ~ 2.5 arc- 
sec centered on the nucleus. The background was obtained 
from a circular annulus with an outer radius ~ 25 arcsec 
centered on the nucleus, but also excluding circular re- 
gions of radius ~ 2.5 arcsec around X-2 and X-3. The 
background was then rescaled to be appropriate for the 
source extraction cell. 

4.1. Temporal Variability 

We find marginal evidence for variability from the nucleus 
of NGC 3226. In Fig. 2a we show the observed count rate 
in 0.8-6 keV band as a function of time using 2 x 10 3 s 
bins. There are 21 such bins, each containing ^20 counts 
on average, thus x 2 -statistics can be used (but with some 
caution). The hypothesis of a constant count rate is re- 
jected at ~ 94% confidence (x 2 = 1-53 for 20 degrees of 
freedom, hereafter d.o.f.). Using 4 x 10 3 s bins, a constant 
count rate is rejected at ~ 97% confidence (xt — 1-96 for 
10 d.o.f.). 

Given the data from NGC 3226 are in the "grey-area" 
where the appropriateness of using % 2 -statistics is debat- 
able (i.e. < 20 bins, each with a relatively small number 
of counts) we have also performed a temporal analysis as- 
suming Poissonian statistics. This analysis also indicates 
evidence for variability. For example, in Fig. 2b we show 
the observed (red) number of bins (each of 10 3 s), contain- 
ing 0, 1, 2, .... counts compared to the predicted Poisson 
distribution (blue) derived from the mean count rate of 
~ 0.01 ct s _1 . For the example shown, we find 6 of the 42 
bins outside the the 95% confidence limits (dashed lines 
in Fig. 2b) compared to the 2.1 predicted from a Poisson 
distribution. We have repeated this analysis 20 times, on 
each occasion offsetting the start time of each 10 3 s bin by 
1/20 of a bin width. We find the observed number of bins 
exceeds the predicted number on 13 of the 20 occasions. 
Similarly we have repeated the analysis using bins of dif- 
ferent temporal sizes, and find cases where number of bins 
observed outside the 95% confidence limits exceeds that 
expected for bins in the ^fewx 10 2 -fewx 10 3 s. Neverthe- 



less, given the current data are in the Poissonian regime, 
again we consider the evidence only marginal. 

4.2. Spectral Analysis 

Spectra were extracted from the (0 th -order) cell described 
above using the psextract script in CIA0. The resultant 
Pulse-Invariant (PI) channels at energies > 0.5 keV were 
then grouped so as to contain a minimum of 20 counts per 
bin, and hence allowing a x 2 minimization technique to be 
employed. This resulted in 24 spectral bins which were an- 
alyzed within the XSPEC (v 11.0.1) package. For the Galac- 
tic absorption along the line-of-sight we assume an effec- 
tive hydrogen column density of Nffj = 2.08 x 10 20 cm~ 2 
(as measured towards NGC 3227; Murphy et al 1996). 

We find an absorbed power law provides an adequate fit 
to the data (% 2 = 20.1 for 21 degrees of freedom). The 
best-fitting parameters and 68% confidence ranges arc a 
photon index T = 1.94±g;||, and intrinsic absorption (in 

addition to N^\) of N H = (4.8±^) x 10 21 cm -2 . The 
implied luminosity (corrected for absorption) in the 0.4— 
10 keV band is L(0.4-10 keV) ~ 3.5±{i x 10 io h^ erg s~\ 
and that in the 2-10 keV band £(2-10 keV) ~ l-8±g;| x 
10 4U /i 

75 er g s 1 - The data and this model are shown in 
Fig. 3a. Unfortunately any Fe Ka emission cannot be well 
constrained by the data (equivalent widths EW < 0.9 keV 
and < 1.5 keV for narrow and broad lines at 6.4 keV 
respectively). An absorbed thermal bremsstrahlung also 
provides an adequate fit to the data (\ 2 — 22.2 for 21 
degrees of freedom), with kT = 6. 51^' i keV, and with a 

luminosity and Njj (in addition to Nffj) consistent with 
those found above. We shall discuss these two models fur- 
ther in §4.3. 

The nucleus of NGC 3226 is sufficiently bright to also 
be detected in the dispersed spectrum. Source and back- 
ground spectra were extracted using the procedure de- 
scribed in George et al. (2001). The source exceeds the 
mean background level in the ±l st -orders in the ranges 
0.18-1.7 nm (MEG) and 0.15-1.0 nm (HEG). The spectral 
resolution of the MEG and HEG can be well approximated 
by Gaussians with ("ler") widths ~ 1 pm and ~ 0.5 pm 
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Fig. 3. — The upper panels show the spectra from the nucleus of NGC 3226, along with the absorbed, power-law model discussed in the 
text. Bins which contain more than 1 photon at > 97.5% confidence are shown in red. The lower panels we show the effective area of the 
instruments, a) The O th -order spectrum grouped such that each bin contains > 20 counts, b) The MEG l st -order spectrum derived using 
only data < 1.6 nm (— l st -order) and < 1.8 nm (+l st -order) to maximize the signal-to-noise ratio. For clarity the data are shown using 23 pm 
wide bins (ie. 10 X the FWHM of the MEG spectral resolution). No statistically significant features are present in the spectrum made at the 
full MEG resolution. 



respectively (1 pm = 10 12 m; and also corresponds to 
-1 ACIS pixel 11 ). Both the MEG and HEG spectra 
have a very low signal-to- noise ratio (< 0.5 ct pm' 1 and 
< 0.3 ct pm -1 respectively). We have constructed spectra 
using bins between 1 and 4 times the resolution of each 
grating arm. However we find no cases where the number 
of counts per bin exceeds 7 (~ 99% confidence limit for 1 
count per bin), and hence no evidence that the spectrum 
contains any intense emission lines. Using larger bin sizes, 
we find both the MEG and HEG spectra are in good agree- 
ment with that derived from the th -order (e.g. Fig.3b). 
Thus we find no requirement for additional spectral com- 
plexity from either the th -order or dispersed spectra. 

Our spectral models predict ROSATHKI count rates (~ 
4x 10" 3 ct s" 1 ) consistent with those seen during an obser- 
vation in 1995 Aug (Roberts & Warwick 2000). However 
the predicted ROSAT PSPC count rates (~ 10" 2 ct s" 1 ) 
are slightly lower than that reported during an observa- 
tion in 1993 May (~4 x 10~ 2 ct s _1 from within a circu- 
lar region of radius ~ ^5h^ kpc is reported by Sansom 
et al. 2000). This may be the result of variability on 
long timescales, and/or an additional soft X-ray compo- 
nent from the AGN. However at least part of the higher 
count rate seen by the PSPC may be due to other sources 
of (perhaps time-variable) X-ray emission associated with 
the galaxy. 

4.3. The Case for a LLAGN at the Center of NGC 3226 

The improved spatial resolution of CXO has enabled us 
to show that the X-ray emission from NGC 3226 is dom- 
inated by a point-like (tentatively variable) source within 
~ 0.2/iyg 1 kpc of the nucleus. The X-ray spectrum is 
measured up to ~10 keV and is consistent with a power 
law with a photon index 1.7 < T < 2.2, or thermal 
bremsstrahlung emission with 4 < kT < 10 keV. The ra- 
tio of observed flux at 2 keV (~ 4 x 10~ 2 [iJy) to that 
at 5 GHz from the flat-spectrum, compact radio source 
(Falcke et al. 2000) gives a radio-to-X-ray spectral index 
a rx = -log(/(2 keV)//(5 GHz))/7.986 = 0.62. Such a 

11 Due to the orientation of the grating arms and the dithering of the 
and wavelength is complex and time-variable. 



value is consistent with that derived from several similar 
objects, including M 81 (NGC 3031) and NGC 4579 (both 
with a rx =~ 0.6; e.g. Ho 1999, and references therein). 
M 81 and NGC 4579 are of note since both are commonly 
believed to be accretion-powered objects. Indeed Fe Ka 
emission has been detected from both objects, which is 
broad in the case of M 81 (Ishisaki et al. 1996) and vari- 
able in the case of NGC 4579 (Terashima et al. 2000), 
highly suggestive of an accretion disk. 

In common with many other LLAGN, optical images 
reveal large quantities of dust in filamentary structures 
within the inner regions of NGC 3226 (e.g. Rest et al. 
2001). The X-ray presented here also indicate intrin- 
sic absorption. The equivalent hydrogen column den- 
sity reported in §4.2 (3-7 x 10 21 cm" 2 ) corresponds to 
2 < Ay < 4 assuming a Galactic dust-to-gas ratio. Using 
the spectral energy distributions of M 81 and NGC 4579 
as templates, both of which have an optical-to- X-ray spec- 
tral index a ox = -log(/(2 keV)//(250 nm))/2.605 ~ 1.0 
(after correcting for the better-known reddening in these 
galaxies), we predict the nucleus has my ~ tur > 21. The 
nucleus is therefore predicted to be swamped by the stellar 
emission (Rest et al. find ~ 14 arcsec -2 ) in the central 
regions, consistent with observations. 

These characteristics strongly suggest NGC 3226 hosts 
a central AGN. Such an hypothesis is also consistent with 
the broad component to the Ha emission line suggested 
by Ho et al. (1997b) based on modeling of a fairly noisy 
spectrum of the Ha region (where Ha is blended with [N 
II] and [S n] ; see their Fig. 10a). Terashima, Ho & 
Ptak (2000) have reported a strong correlation between 
L(2-10 keV) and the luminosity of Ha for a number of 
galaxies containing type 1 LINERs. This correlation has 
recently been confirmed and extended by Ho et al. (2001). 
This correlation supports the hypothesis that the domi- 
nant ionization source of the type 1 LINERs is the photon 
field of a central LLAGN. NGC 3226 is consistent with 
this correlation, adding further support to this hypothe- 
sis. The fact that the observed broad component of Ha 

spacecraft during the observation the relationship between ACIS pixel 
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constitutes a lower faction (~ 60%) of the total (broad 
plus narrow) line emission compared to a more "normal" 
Seyfert 1 galaxy (eg. <~ 97% in NGC 3516) suggests differ- 
ential (probably patchy) absorption towards the narrow- 
and broad-line regions. 

4.3.1. Constraints on the Accretion Process 

The optical luminosity of the entire galaxy NGC 3226 
(M Bt = -19.4; L B /L Q ~ 8 x 10 9 ) and a central, line-of- 
sight velocity dispersion (~ 180 km s _1 ; Simien & Prug- 
niel 1998) are consistent with the relation found for a va- 
riety of galaxies harboring black holes (eg Gebhardt et al. 
2000a; and references therein). The implied central black 
hole mass is ~ 10 8 M©, and hence the implied accretion 
rate is m ~ 4 x 10 _6 /i 75 1 / b oix M Edd , where / bo ix is the 
ratio of bolometric to 0.4-10 keV luminosity (which is un- 
likely to much exceed 10), and M Edd is the Eddington ac- 
cretion rate. Some words of caution are necessary. First, 
we note that the line-of-sight velocity dispersion quoted 
above is derived from ground-based observations - as yet 
there have been no HST observations capable of resolving 
the central kinematics of NGC 3226. Second, there have 
recently been a number of questions raised regarding the 
exact relation between the luminosity of the bulge and the 
mass of the black hole estimated by various means in dif- 
ferent samples (eg. see Gebhardt et al. 2000b; Merritt & 
Ferrarese 2001, and references therein). Nevertheless, if a 
black hole indeed exists at the center of NGC 3226, the 
accretion rate is surely rii< 10~ 2 -MEdd- 

In systems with such low rates, the mass infall is gen- 
erally thought to occur via a "standard" (geometrically- 
thin, viscosity-dominated) accretion disk at large radii, but 
switch to a radiatively less-efficient advection-dominated 
accretion flow (ADAF) in the inner regions (e.g. Narayan, 
Mahadevan & Quataert 1998, and references therein). The 
transition radius at which the flow becomes advection- 
dominated depends on a number of model parameters. 
Estimates range from ~ 10-10 2 Schwarzschild radii for 
several LLAGN (including M 81 and NGC 4579; e.g. see 
Quataert et al. 1999) to ~ 10 3 Schwarzschild radii for 
giant elliptical galaxies within clusters (Di Matteo et al. 
2000). The X-ray emission from ADAFs is believed to 
be dominated by thermal bremsstrahlung (with kT ~few 
keV), and not expected to exhibit significantly variabil- 
ity on timescales shorter than the local dynamical time 
(~ 5 x 10 4 s at 10 Schwarzschild radii from a 10 8 M© 
black hole - see also Ptak et al. 1998). In §4.2 we found 
a thermal bremsstrahlung model to be consistent with the 
observed spectrum of the nucleus. However, if the vari- 
ability on timescales <fewxl0 3 s suggested in §4.1 is con- 
firmed, time-dependent models will be required. 

5. THE OFF-NUCLEAR SOURCES CXO J102334. 1 + 195347 

and exo J102326. 7+195407 

There are too few counts from cither CXO J102334. 1+195347 
AND CXO J102326.7+195407 to enable detailed analysis. 
We find no evidence for statistically significant variability 
on any timescale. Nevertheless some crude insight can be 
gained from consideration of counts detected in different 
energy bands. For CXO J102334.1+195347, the ratio of 
source counts (0.3-2 keV)/(2-10 keV) = 0.84 + 0.40. For 
a power law spectrum, such a ratio can be obtained with 



0-9 < T < 1.8 if absorption by only Nf^\ is assumed, 
or with r = 2 and additional absorption of Nh — 10 21 ~ 
10 22 cm" 2 . In the case of CXO J102326. 7+195407, we 
find the ratio of counts (0.3-2 kcV)/(2-10 keV) < 0.2, 
indicative of an even harder spectrum. Assuming absorp- 
tion by only Nf^\, a power law with Y < 0.5 is required. 
Assuming a power law with T = 2, additional absorption 
with Nh > 2 x 10 22 cm -2 is required. 

The proximity of both sources to the nucleus (within 
a projected distance ~ 1.3/ifg 1 kpc) strongly suggests 
both are X-ray sources within NGC 3226. The observed 
flux from each source is ~fewxl0~ 14 erg cm -2 s _1 in the 
2-10 keV band. The density of background sources is 
< 10 2 deg -2 in this band at this flux level (e.g. Mushotzky 
et al 2000; Tozzi et al. 2001; and references therein). Thus 
the probability of a background source this bright within 
15 arcsec of the nucleus of NGC 3226 is < 10~ 4 . 

Assuming both sources arc indeed within NGC 3226, 
the spectral forms described above imply absorption- 
corrected luminosities in the range ^fewxlO 38 
fewx lO 39 ^^ 1 erg s _1 . Assuming all emission from each 
source is indeed due to a single object, either or both 
CXO J102334.1+195347 and CXO J102326. 7+195407 
therefore have luminosities exceeding the Eddington limit 
for neutron stars (~ 3 x 10 38 erg s _1 ), and hence must 
be >fewM© black-holes, very young supernovae or micro- 
quasars. Such super luminal sources are now being found 
regularly by CXO (eg. Griffiths et al. 2000; Fabbiano, 
Zezas & Murray 2001). Variability and detailed X-ray 
spectroscopic studies are the most likely means the vari- 
ous source types can be distinguished. 



6. CONCLUSIONS AND FUTURE PROSPECTS 

The results presented here illustrate the power of X-ray 
observations to detect and characterize the emission from 
LLAGN: a class of source which can usually only other- 
wise be studied at radio wavelengths and possibly with 
high signal-to-noise, high-resolution optical emission-line 
studies. Indeed in many cases, as for NGC 3226, the 
nucleus is expected to be invisible to HST in the opti- 
cal band. We anticipate that many such objects will be 
detected serendipitously in the fields of future CXO and 
XMM-Newton observations. These data will expand the 
regions of parameter space accessible to study, surely re- 
sulting in a better understanding of the AGN/black-hole 
phenomenon. 
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